ON THE NATURE AND THE ORDER OF THE DECONFINING 

TRANSITION IN QCD. 
In honour of Yu. A. Simonov on his seventyth birthday. 



Adriano Di Giacomo^ 

^ Dipartimento di Fisica dell'Universita di Pisa , and INFN Sezione di Pisa 
i-j.". l.go Pontecorvo 2, 1-56127, Pisa, Italy 

o . 

■ The determination of the parameters of the deconfining transition in Nf = 2 

bJQ. QCD is discussed , and its relevance to the understanding of the mechanism of color 

I confinement. 
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1. INTRODUCTION. 



■ Understanding the mechanism by which QCD confines colored particles is one of the most 

00 ; 

O . fundamental and challenging problems in field theory [1] [2]. 

O I Experiments aimed to detect the existence of a deconfining transition at high temperature 

ci . by colliding heavy ions have not yet produced a definite answer [3]. Up to now the decon- 
Qh! fining transition has been only observed in numerical simulations of QCD on a lattice[4]. 

' Identifying a signal of deconfinement is far from trivial in experiments, and is not easy in 
- ^ , lattice simulations either. To be rigorous confinement means absence of colored particles 
' in asymptotic states. Of course one cannot check confinement by looking at all asymptotic 
states. In the quenched approximation, in which dynamical quark loops are neglected, one 
looks at the static potential acting between a heavy Q — Q pair. The criterion for confinement 
is then the existence of a string tension, which means a linear behaviour at large distances 

V{r) = ar (1) 

The static potential is measured on the lattice either in terms of the Wilson loop (parallel 
transport along a closed path), or equivalently in terms of the correlator of Polyakov lines 
(parallel transports along the time axis). For a Wilson loop of extension t in the time 
direction and r in space 

W{r,t) oc exp[-tV{r)] (2) 
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Confinement eq(l) implies the so-called area law 

W{r,t) (xex.p[-art] (3) 

r, t — > cx) means that they are large compared to the correlation length. The correlator of 
two Polyakov lines can be written, by use of cluster property, 

(L(r)Lt(O)) ~ Cexp[-ar/T] + \{L)f (4) 

) — >oo 

On the other hand 

V(r) ~ -T\n{L{r)L\0)) (5) 

It is found by numerical simulations that a critical temperature Tc exists such that for T < Tc 
(L) =0. From eq's(4) and (5) confinement eq(l) then follows. For T > Tc (L) ^ and 

V{r) — >■ const (6) 

r,t— >oo 

For quenched SU{3) T^j -sfla) ^ .65 [5]. A deconfining phase transition exists at Tc ~ 
270Mef and the Polyakov line (L) is the order parameter. The symmetry involved is Z3 
which is broken at T > Tc, and the transition is order-disorder. Some problem can arise in 
the continuum limit with the definition [6]of the order parameter (L) , but the main idea 
looks sound. 

The situation is less clear in the more realistic case of full QCD, including dynamical 
quarks. There Z3 is not a symmetry anyhow, being broken by the coupling to the quarks. 
Moreover the potential at large distances is not expected to grow with r, due to the conver- 
sion of the potential energy into light q — q pairs (string breaking) [7]. One needs then an 
alternative criterion for confinement, which , however, is not known. 

A phase transition is expected anyhow at low quark masses from the low temperature 
phase in which chiral symmetry is spontaneously broken, to a phase in which it is restored: 
the order parameter for this transition is the chiral condensate (^V"). This transition is 
indeed observed on the lattice at ~ 170Mev. A priori restoration of chiral symmetry is 
not the same as deconfinement, even if it is physically understandable , e.g. by thinking 
of a bag model, that confinement can imply chiral symmetry breaking. In the absence of 
a criterion for confinement the question if chiral and deconfining transition coincide can 
not even be asked. A third symmetry exists, the axial U(l) , which is broken at low T 
as an effect of the anomaly, and is also expected to be restored at high temperature. In 
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principle the corresponding transition temperature could be different from that of the other 
two transitions. 

An analysis can be done of the chiral transition based on symmetry arguments and on 
renormalization group techniques [8]. If the assumption is made that the relevant degrees 
of freedom at the transition are the pseudoscalar particles an effective lagrangean can be 
written on the basis of symmetry and scale invariance, which describes the density of free 
energy around the chiral point. For Nf > 3 no infrared stable fixed point exists , so that 
the transition is expected to be first order. Nf — 2 is a special case, and the transition can 
either be weak first order or second order . If it is first order , it is weak first order also at 
rriq 7^ 0; if it is second order the transition at nig 7^ is a crossover [8]. What is the case 
can be investigated by numerical simulations on the lattice, even if up to now [9] [10] [11] the 
results have been rather elusive. 

This issue is very important for the understanding of the deconfining phase transition. On 
the basis of the quenched case mentioned above the deconfining transition should be order- 
disorder, and a genuine order parameter should exist labeling the two phases. Nc ^ 00 
arguments suggest indeed that the symmetry involved as well as the mechanism of confine- 
ment should be the same for quenched and unquenched and independent. If for Nf = 2 
the chiral transition is second order, then the transition for rUg 7^ is a crossover, and the 
deconfining transition is not order-disorder. A first order chiral transition would instead be 
consistent with order-disorder, and possibly be such up to niq — 00 , which is the quenched 
case. 

The order of the transition can be investigated by looking at the volume dependence of the 
specific heat in numerical simulations, by a technique known as finite size scaling [see sect 
2 below]. The free energy density around the phase transition is determined by symmetry 
arguments [12] up to unknown numerical coefficients, in terms of the order parameter. The 
(pseudo) critical indices determined by looking at the susceptibility of the order parameter 
must agree with the determination made by looking at the specific heat . Such an agreement 
is needed to legitimate the choice of the order paramenter. 

An order parameter {^) has been developed and tested [13] [14] [15], based on the working 
hypothesis that the mechanism of confinement is dual superconductivity of the vacuum[2]. 
(/x) is the vacuum expectation value of an operator carrying magnetic charge, and , unless 
the Polyakov loop or the chiral condensate is well defined independent of Nc,Nf even in the 
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continuum limit. 

The critical indices can be investigated by measuring susceptibilities involving /i , in 
particular the quantity p = (i/(ir[ln((/i))] , with r the reduced temperature 

T = 1 - T/Te (7) 

If they agree with the determination made by use of the specific heat, an additional le- 
gitimation results for the order parameter and for the mechanism of confinement by dual 
superconductivity. Confinement will then be defined in terms of an appropriate symmetry, 
and the question whether the deconfining transition and the chiral transitions occurr at the 
same temperature becomes meaningful. 

In sect 2 we shall discuss new results on the above issues[16],[17] . 

In sects we draw some conclusion. 



2. LATTICE RESULTS AND FINITE SIZE SCALING ANALYSIS. 

The theory of finite size scaling for higher order and weak first order phase transitions 
is based on renormalization group arguments [18] [19] , which are expected to hold when, 
by approaching the transition, the correlation lengths become much larger than the lattice 
spacing. It allows to extrapolate to infinite volume results obtained at finite volumes. The 
extrapolation depends on the (pseudo) critical indices , which can then be determined from 
the volume dependence, and with them the order and the universality class of the transition. 

The relevant quantities , which are related to derivatives of the free energy, and hence to 
the order of the transition, are the susceptibilities . In the following we shall consider the 
susceptibility of the chiral order parameter 

T 

the specific heat 

1 d'^ 

and 

p=|:ln(/i) (10) 

All these susceptibilities will depend on the parameters which describe the system, i.e. 
P — ^ with g the gauge coupling, and ma , the bare quark mass in units of the inverse 
lattice spacing. 
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The simulation is made on a lattice NtN^^Nt <^ Ns, of extension Nt spacings along the 
time axis and Ns along the three spatial directions, with periodic boundary conditions (be) 
in time for bosons, antiperiodic be for fermions. 

The temperature is the inverse of the time extension, or 

(11) 

Nta{f3,ma) ^ ' 

a(/?, ma) is the lattice spacing in physical units. 

All the susceptibilities as function of the temperature will have peaks at Tc, which diverge 
as the volume V ^ with a power depending on the order of the transition , specifically on 
the critical indices. In the usual notation of statistical mechanics, denoting by L the spatial 
size of the system 

X = L^/''<^^{TL^,mLy^) (12) 
Cv = c^ + L"/'^$e(TL^,mL^'') (13) 
p = L^^^{TL^,mLy'^) (14) 



Since T = „ ,\ — r 

Nta(P,ma) 



a{Pc,0) ^^^^ 



a{(3, ma) 

In the vicinity of the chiral point one can expand in powers of /3c — /3 and of ma, obtaining 



T (X /3c - /3 + Kma (16) 

with 

rl]r\ n rl]r\ n 

(17) 



^ d\na ^dlna 



d{ma) d/3 

In the quenched case the last term in eq(16) is absent. 

When the volume of the system goes large at tL~ fixed, the susceptibilities will tend to a 
finite limit certainly if the transition is second order, and the free energy is continuous and 
finite in the neighbourhood of the chiral point. This implies that the powers of L in front of 
the susceptibility eq's (12), (13), (14) have to be eliminated by the dependence on mW^ , and 
this means for the three susceptibilities 

X oc m~~^F^{TL^) (18) 

n — 1 

Cv — Cy OC m •'yhFcirL") (19) 
p (X m ''yhFp{TL-^) (20) 
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If the transition is weak first order the system will behave in the same way for intermediate 
values of L, smaller or equal than the critical correlation length, apart for the different value 
of the critical indices. 

Eq's(18,(19),(20) imply that the maxima of the susceptibilities with respect to r lie on 
the line of the plane ma) of equation 

tL^ = C (21) 

with C a constant, or, by use of eq(16) 

C 

Pc-P + Kma r = (22) 

This relation can be tested against the lattice data , to determine the value of v . We 
will discuss the result in the following . 

Eqs (18). (19), (20) also imply that height of the the maximum of the susceptibilities is 
volume independent, whatever the order of the transition, in the limits of the regime in 
which these equations are obeyed when the transition is first order. 

For second order transition the result is valid for any volume and fits the idea that the 
transition line in the plane (/?, ma) is a crossover. 

For a first order transition ,if we cross the transition line at fixed ma, we expect a be- 
haviour typical of a first order transition, namely a growth of the peaks of the susceptibilities 
proportional to the volume. Indeed for a first order transition a = 7 = l,z/ = 1/d = 1/3. 
This behaviour will be visible when L becomes larger than the critical correlation length, 
and together with it the typical two-peak distribution of the internal energy will appear. 

Moreover, if we keep mL^** constant, eqs(12),(13),(14)) require that the maximum of the 
peaks is again proportional to a well defined power of L depending on the nature of the 
transition. 

The details of the simulations will be reported elsewhere [16] [17] Here we only summarize 
the results. 

As a strategy we assume the critical indices expected for second order 0(4) phase tran- 
sition and in alternative those of a first order transition and we investigate wether data are 
consistent with either of them. A measure of the agreement is the value of x^/dof. 

In table 1 we recall values of the critical indices for the two cases. 

We investigate 
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1) The position of the peak as a function of mass and volume , which scales according 
to Eq(22). The ^jdoj is ~ 1 for the choice first order, and typically ~ 10 for 0(4) 
second order. 

2) The peak values of the susceptibilities. They all occurr on the same line, within 
errors. The height of the peak is volume independent at fixed ma for moderate values 
of the volume and proportional to an appropriate power of (ma) Eqs(18)-(20), which 
depends on the nature of the transition. Here again the y^jdoj is compatible with 1 
for a first order transition , and much worse , typically [ /dof ~ 10] for 0(4) second 
order. 

3) For low values of m (e.g. ma = .0135 ) as the volume increases further a growth of 
the height of the peak is observed and some sign of bistability in the time histories of 
the energy density. An increase of the heigth of the peak by a factor ^ 2 is observed 
going from 16^ to 32^ . This means that there is a transition, and not a crossover. At 
higher values of ma the transition looks weaker and we were not yet able to reach a 
large enough volume to see a bistability . 

4) The behaviour with respect to volume at fixed mL^*" is consistent with first order 
transition , and disfavours second order : typically jdoj ~ 1 for first order, jdoj ~ 
10 for second order 0(4). 

Previous investigations of the same system were made on rather small volumes[10][9] [11] or 
in a less systematic way, and the results where admittedly inconclusive, although a slight psi- 
cological preference appeared for a second order transition . Numerically our determinations 
are consistent with the previous ones when done at the same values of the parameters. The 
details of the simulations and the comparison with previous work will be presented elsewhere 
[16][17]. 

3. DISCUSSION. 

We have obtained substantial preliminary evidence that the chiral transition for Nf — 2 
QCD is first order. This makes the deconfining transition consistent with an order disorder 
transition. The disorder parameter {ji) detecting dual superconductivity of the vacuum 
provides a determination of the critical indices consistent with the one based on the specific 
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heat, and a transition line consistent with the other susceptibilities. It can then be used as an 
order parameter for confinement. The deconfining transition and the chiral transition occurr 
at the same T. A careful numerical study of the U(l) axial anomaly across the deconfining 
transition is on the way and will provide a check of the assumption made in [8] that the 
relevant degrees of freedom at the chiral transition are the pseudoscalar goldstone particles. 

I wish to thank all my collaborators , M. D'Elia, B.Lucini, G. Paffuti, C. Pica who 
contributed substantially to the results presented. Many ideas and motivations also come 
from discussions and collaboration with Yu. A. Simonov , to whom this paper is dedicated 
with deep friendship in the occasion of his 70th birthday. 
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Table 1. 

The (pseudo) critical indices for weak 1st order and second order 0(4) phase transition. 





a 


7 


u 


Vh 


1st order 
0(4) 2nd order 


1 

-.25(1) 


1 

1.47(1) 


1/3 
.75(1) 


3 

2.49(1) 



